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Abstract. We present the inclusive invariant-mass, transverse momentum and rapidity distributions of 
dielectrons (e + e~ pairs) in p+p interactions at 3.5 GeV beam kinetic energy. In the vector-meson mass 
region, a distinct peak corresponding to direct uj decays is reconstructed with 2% mass resolution. The 
data is compared to predictions from three model calculations. Due to the large acceptance of the HADES 
apparatus for e + e~ invariant masses above 0.2 GeV/c 2 and for transverse pair momenta pt < 1 GeV/c, 
acceptance corrections are to a large extent model independent. This allows us to extract from dielectron 
data for the first time at this energy the inclusive production cross sections for light vector mesons. Inclusive 
production cross sections for tt° and n mesons are also reported. The obtained results will serve as an 
important reference for the study of vector meson production in proton-nucleus and heavy-ion collisions. 
Furthermore, using this data, an improved value for the upper bound of the branching ratio for direct n 
decays into the electron-positron channel is obtained. 
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1 Introduction 

The High Acceptance Di-Electron Spectrometer (HADES) 
[1 is operated at the GSI Helmholtzzentrum fur Schweri- 
onenforschung in Darmstadt, Germany. One of the main 
physics goals of HADES is to investigate spectral modifi- 
cations of light vector mesons in strongly interacting mat- 
ter. The question is how the low-energy QCD spectrum, 
which is experimentally known in the vacuum in terms of 
hadron spectra, will change when this vacuum is heated 
and filled up with color charges. Spectral modifications of 
hadrons (encoded, e.g., by changes of their masses and de- 
cay widths) in hot and/or dense matter are often discussed 
in the context of the restoration of the broken chiral sym- 
metry. Detailed investigations, however, reveal that the 
link between hadron properties and QCD symmetries is 
not as direct as originally envisaged [2[3]. 

In order to search for in-medium effects, results on 
electron-positron invariant mass spectra from proton in- 
duced reactions on nuclei and from heavy-ion collisions 
should be compared systematically and complemented with 
insights gained from photo-induced reactions. The inter- 
pretation of nuclear data requires firm knowledge of the 
corresponding data from proton-proton (p+p) collisions. 
These investigations are also important for transport model 
calculations, as the dilepton spectra from elementary in- 
teractions serve as input to these codes. Given this motiva- 
tion, the HADES collaboration has set up an experimental 
program to measure dilepton spectra in elementary colli- 
sions. First results were discussed in [4]. Here we present 
the inclusive dielectron spectra measured in p+p collisions 
at 3.5 GeV kinetic beam energy. 

At this energy the dominant mechanisms for hadron 
production are still controversially discussed and inclu- 
sive production cross sections are not yet measured. In 
fact, various production scenarios are assumed in differ- 
ent transport codes: While in HSEQ [5! and GiBUlQ [5] 
hadrons are produced at this beam energy through string 
fragmentation \7\ , in UrQMEP] [5] the decays of nucleon 
resonances are the sources of final state particles, in par- 
ticular of mesons. 

Another uncertainty concerns the Dalitz decays of baryon 
resonances (R — > N e + e~). Here two different aspects are 
a matter of discussion. First, the R Dalitz decay process, 
i.e. its decay into a nucleon and a massive photon with the 
subsequent decay of the latter into a dielectron depends 
on the electromagnetic structure of the N-R transition ver- 
tex. In case of the Z\(1232) resonance this transition vertex 
depends on three independent helicity amplitudes corre- 
sponding to three helicity states of the massive photon 
and two of the nucleon. Equivalently, one can describe 
this vertex by three independent transition form factors 
built up from the helicity amplitudes. In this decay pro- 
cess the squared four-momentum of the virtual photon 
equals the squared invariant mass of the lepton pair and 
therefore is a positive quantity. Such a photon transfers 



energy and is said to be a time-like photorj^J There are 
also other processes (e.g. pion electro-production) where 
the virtual photon transfers momentum, but not energy, 
and is then referred to as space-like. While in case of 
space-like photons the above mentioned transition form 
factors have been measured in quite a wide range of four- 
momentum transfer, for time-like photons their mass de- 
pendence is not settled yet. The mass dependence of the 
form factors is usually modeled within the Vector Me- 
son Dominance (VMD) model. In this picture the virtual 
photon couples to the nucleon through intermediate vec- 
tor meson (p, w, </>) states. However, it is known that the 
standard VMD model of Sakurai [SlUU] significantly over- 
estimates the radiative (R — > N7) branching ratios once 
the coupling constant is extracted from the corresponding 
mesonic decays. The modified VMD model of Kroll et al. 
[TT] allows to fix the pN and 7N coupling constants inde- 
pendently. The asymptotic behavior of the transition form 
factors predicted by both models disagrees however with 
the outcome of quark counting rules [12]. On the other 
hand, the authors of [T5] demonstrated in their frame- 
work of extended VMD that the incorporation of higher 
vector meson states resolves the problem between photon 
and p meson branchings of the nucleon resonances. Yet 
another VMD model variant for the nucleon form factors 
was proposed by Iachello [13]. It describes simultaneously 
the nucleon space-like and time-like form factors as well 
as space- like N-Z\ transition form factors [T4"lll5| . 

A second aspect of resonance decays to consider is the 
parametrization of the mass dependent resonance width. 
Various prescriptions are used in model calculations (for 
details see section 3 and [16]) which differ at high reso- 
nance masses. As a consequence, the resulting dielectron 
yield from the resonance Dalitz decays has a large uncer- 
tainty, as shown for the Zi(1232) in [16]. 

Our precision data offer hence a unique possibility to 
address the above-mentioned problems. 

This paper is organized as follows. After giving detailed 
information about the collected data and analysis chain 
in section [2] we report in section |on a comparison of 
the data with model calculations. The sensitivity of the 
data to the N-Zi electromagnetic transition vertex is also 
discussed in this section. The extraction of 7r°, 77, Z\(1232), 
p and ui cross sections is discussed in section [4] The data 
allow to improve the upper bound of the direct rj — > e + e~ 
decay, as discussed in section [5] 



2 The data 

In the experiment (see pQ for a detailed description of 
HADES), a proton beam of 10 7 particles/s with a kinetic 
energy of 3.5 GeV was incident on a 4.4 cm long liquid hy- 
drogen target [17]. The data readout was started upon a 
first-level trigger (LVL1) decision. Depending on the re- 
action channel of interest, two different settings of the 



Hadron-String Dynamics 

2 Giessen Boltzmann-Uehling-Uhlenbeck project 

3 Ultrarelativistic Quantum Molecular Dynamics 



We use the following convention of the squared 4-momenta: 
P 2 = P fl (E,p)Pii(E, —p) = E 2 — \p\ 2 , i.e., positive p 2 is called 
time-like, and negative p 2 is correspondingly space-like. 
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LVL1 trigger were required: (i) a charged-particle mul- 
tiplicity MULT > 3 to enhance inclusive dielectron pro- 
duction and (ii) MULT > 2 with hits in opposite sectors 
{A<f> = 180° ± 60°) of the time-of-flight detectors to en- 
rich elastic p+p events used for the absolute normaliza- 
tion of the dielectron data. The LVL1 was followed by 
a second-level trigger (LVL2) requesting at least one elec- 
tron candidate recognized in the Ring-Imaging Cherenkov 
Detector (RICH) and time-of-flight /pre-shower detectors 
[JJ . All events with positive LVL2 trigger decision and ev- 
ery third LVL1 event, irrespective of the LVL2 decision, 
were recorded, yielding a total of 1.17xl0 9 events. 




M ee [GeV/c 2 ] 

Fig. 1. (color online) Efficiency corrected inclusive invariant- 
mass distribution of dielectrons inside the geometrical accep- 
tance of HADES of signal (black dots) e + e~ pairs after combi- 
natorial background (blue dots) subtraction for p(3.5GeV)+p 
collisions. In the right upper part of the figure the uj meson 
region is shown on a linear scale. The data are normalized to 
the simultaneously measured p+p elastic events. 

The identified single-electron and single-positron tracks 
were combined into opposite-sign pairs. However, many 
of these pairs represent combinatorial background (CB) 
which is mostly due to uncorrelated pairs from multi-pion 
Dalitz decays and correlated ones from ir° — > 77 accompa- 
nied by photon conversion in the detector material and/or 
from Dalitz decays. The combinatorial background was re- 
duced by using an opening angle cut of a ee > 9° between 
the reconstructed lepton tracks and a condition on quality 
criteria of the track fitting algorithm [1] . In addition a mo- 
mentum cut of 0.08 < p[GeV/c] < 2 for each lepton was 
applied. The combinatorial background was formed from 
the sum of the reconstructed like-sign invariant-mass dis- 



tributions, dN ++ /dM ee and dN — /dM ee . The like-sign 
pairs were subjected to the same selection criteria as the 
unlike-sign ones. Both, the unlike sign invariant-mass dis- 
tribution and the CB, were corrected for the detector and 
reconstruction inefficiencies on a pair-by-pair basis, de- 
fined as the product of single-lepton efficiencies deduced 
from dedicated Monte-Carlo events embedded into real 
events [JJ[T8, 19 • The geometrical pair acceptance (the ac- 
ceptance matrix) of the HADES detector was obtained in 
a similar way from single-lepton acceptances defined as a 
function of the lepton momentum, polar and azimuthal 
emission angles pQ. 

The final signal pair distribution inside the geomet- 
rical acceptance of the HADES, shown in Fig. [l] is the 
result of a subtraction of the CB from the unlike-sign 
invariant-mass distribution. Both spectra are normalized 
to N"f c /a"i c , where N^ cc and a®f c denote the measured 
yield of the p+p elastic scattering and the differential elas- 
tic cross section [20] inside the acceptance of HADES, re- 
spectively. 

The low-mass region of the spectrum in Fig. IT] is dom- 
inated by Dalitz decays of neutral mesons (tt° , 77, uj), as 
well as by the Dalitz decay of the Zi(1232) resonance (see 
below). The evident peak around the pole mass of the 
u) meson corresponds to its direct decay into e + e~ pairs. 
However, this mass range contains also pairs stemming 
from the direct decays of the p meson. In total, 6.1 xlO 4 
signal pairs, 5.4xl0 4 of them in the region below 0.15 
GeV/c 2 , were reconstructed. The number of pairs in the 
mass range between 0.71 GeV/c 2 and 0.81 GeV/c 2 , which 
corresponds to the ±3cr interval around the reconstructed 
uj peak, amounts to 260. In the inset of Fig. [TJ the uj peak 
is shown on a linear scale. In order to estimate the mass 
resolution, this peak was fitted with a Gaussian distribu- 
tion plus a polynomial for the underlying continuum. The 
obtained mass resolution, c/M" le , is about 2%. 



3 Comparison to three models 

3.1 PYTHIA+PLUTO, UrQMD, HSD 

For better understanding of the inclusive hadron produc- 
tion in 3.5 GeV p+p interactions, we compare in this sec- 
tion the experimentally measured distributions to the re- 
sults from the PYTHIA [21], UrQMD [HIH] and HSD 
[23] event generators. 

PYTHIA, as well as HSD at this energy, use a Monte 
Carlo realization of the Lund string fragmentation model, 
where the assumption of a linear confinement potential 
between the quark and antiquark is taken as a starting 
point. Although the latter is usually used to describe the 
multi-particle production in the high-energy regime, it was 
successfully applied, with some additional adjustments, 
called tunes, to reproduce the experimental data at low 
energies. In particular, since PYTHIA is not predictive in 
assigning the spin to the newly created quark and anti- 
quark pairs from two adjacent string break-ups, it has a 
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tunable parameter which can be adjusted in order to get 
vector meson multiplicities in accordance with our mea- 
sured vector meson yields. Recently such tunes have been 
obtained by the Giessen group, in particular for the p+p 
data at 3.5 GeV jMJUS]. Our tuned values are the same, 
apart from two parameters 

A large fraction of the hadrons produced by fragmenta- 
tion are unstable and subsequently decay into final states. 
We do not let the particles decay directly inside PYTHIA, 
but rather decay them using the PLUTO [26] code. The 
input information obtained from PYTHIA consists thus 
of the particle multiplicities and their four-momenta. 

The simulated dielectron spectra (referred to as the 
cocktail) can be expressed as the incoherent sum over var- 
ious sources of dielectrons, such as Dalitz decays of the 
pseudoscalar mesons n° and ry, Dalitz decay of the vector 
meson uj, Dalitz decay of the 4(1232) as well as direct 
vector meson decays V— > c + c~ with V=p, uj. The plain 
brcmsstrahlung contribution is expected to be small at the 
present beam energy [37] . Similar cocktails have also been 
considered in previous attempts to describe the dielectron 
production in pp, pd and pA collisions [2~8,29,30,3fll3"2l 
33,23,34,25 . Schematically, the differential distribution 
of dielectrons with invariant mass M ee can be expressed 
as a superposition of the above mentioned decay channels 



da 
dM el 



r} oi dM„ 



(i) 



This expression means that a parent hadron i is created in 
the p+p collision with cross section <7j and decays subse- 
quently, thus generating the distribution dTi/ 'dM ee . The 
factor l/rf ot is the inverse of the total width of hadron 
i; together with the partial width for the dielectron de- 
cay channel it encodes the branching ratio. Broad res- 
onances, such as the p and 4(1232), are actually gen- 
erated at masses m p and (see (2) below), and the 
decay distribution depends correspondingly on m p> A and 
M ee (cf. equations (20, 21) in [23]). Note that non-strange 
baryon resonances besides the 4(1232) are not included 
in PYTHIA. Moreover, in case of Dalitz decays, the mass 
dependence of the electromagnetic transition form factors 
should be considered. The mass dependences of the elec- 
tromagnetic transition form factors for the tt°, rj and to 
mesons are parametrized in PLUTO [55] according to [33] 
152"] in agreement with recent measurements [561137] , The 
direct decays of vector mesons are treated within the VMD 
model [5B], while the formulas for the pseudoscalar (P = 
7r° or rj) and vector meson Dalitz decays are adopted from 
[351I39]. The Dalitz decay of the 4(1232) resonance is sim- 
ulated using the expression for its differential decay rate 
from [39 . In these calculations, the N — 4 transition ver- 
tex is described by electric, magnetic and Coulomb form 
factors, corresponding to three independent helicity am- 
plitudes [JO] • As mentioned in the introduction, in the 
time-like region the q 2 dependence of the transition form 
factors is not measured yet. Therefore, we make an ap- 



° In our simulation the values PARP(91) = 
PARJ(21)= 0.36 have been used. 



0.44 and 



proximation by fixing the form factors at the photon point 
(real photons) using the measured radiative decay width 
of the 4(1232) (7^^=0.61-0.7 MeV) gU[3J. We fur- 
ther neglect the terms with the electric form factor, as the 
electric transition is much weaker than the magnetic one 

It should be further noted that in PYTHIA the p and 
4(1232) resonances are implemented with constant total 
widths around the resonance pole mass. This treatment is 
not precise enough outside the resonance pole. Therefore, 
following the prescription of [13115] , we generate masses of 
4(1232), p and uj states inside PYTHIA according to the 
relativistic Breit-Wigner distribution 



A(M) = N 



M 2 r tt 



7T(M 2 



(Mr tot {M)f 



(2) 



with mass dependent total width r tot (M) in case of 4(1232) 
and p, and a constant total width at the pole for the nar- 
row uj meson state. This mass dependence of the total 
width in case of the 4(1232) baryon resonance is calcu- 
lated from its dominant decay channel into pion and nu- 
cleon final states with the cutoff parametrization of [44] 



r t A ot {M A ) 



r 



A—tnN 



r, 



M, 



pole 



pole 



Ma \q po ie 



Qpole 



6 2 



(3) 



where 5 = 0.197 GeV, is the actual mass of the 4(1232), 
M po ; e is its pole mass and r voie is its pole widthQ Further- 
more, q and q po ie denote the pion 3-momenta m the rest 
frame of the 4(1232) with mass and M po / e , respec- 
tively. In [T5] the authors investigate the effect of different 
cutoff prescriptions [44,45 , in particular on the dielec- 
tron spectra for high masses. The resulting uncertainties 
are larger than a factor of 3. 

The mass dependence of the p total width is parametrized 
according to [35] . 

The normalization constant in Eq. AT, is chosen 
such that f™* A(M)dM = 1, where max is fixed at 2 
GeV and min is taken to be 2m 7r , 3m„. and m n +niN for 
p, oj and 4(1232) correspondingly. 

In contrast to PYTHIA and HSD, the transport model 
UrQMD [M] uses a resonance (R) excitation mechanism 
for the production of particles via two-nucleon (N) reac- 
tions of the type NN — >• NR, NN -> RR. Resonances with 
masses up to 2.2 GeV for the N* and 1.95 GeV for the 
4 are included }22j . The production matrix elements for 
resonances are obtained from the experimental data on 
7r, rj and p production, when available. The probabilities 
for the resonances to decay into specific channels arc then 
given by the corresponding known branching ratios. Even 
though in UrQMD the excitation of many baryon reso- 
nances are used, only the Dalitz decay of the 4(1232) 
isobar is explicitly included. The e + e _ contribution from 
decays of higher-lying resonances is included via their p 
decay branches (see table 3.4 in [22]). However, this ap- 
proach leads to an overestimation of the e + e~ production 

6 In this paper c=l units are used for the formulas. 
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from the p decays. On the other hand, estimates of the 
e + e _ yield from the Dalitz decays of higher resonances in- 
dicate smaller contribution as compared to the one from 
4(1232) and r, decays for M ee < 0.55 GeV/c 2 gSEZU^]- 
Dalitz decays of higher-lying baryon resonances have also 
been investigated in [T^] . Checking the validity of this ap- 
proach for describing the experimental data is a subject 
of ongoing HADES activities. 

In general, the reaction pp — > c+ c~ X is fully described 
by three independent degrees of freedom (neglecting the 
internal degrees of freedom, like helicity angles of virtual 
photons) which can be selected in a variety of ways. It is 
important that a given event generator describes the ex- 
perimentally measured distributions in all degrees of free- 
dom. Therefore we present here the comparison of pair 
invariant mass, transverse momentum and rapidity distri- 
butions to the corresponding distributions from simulated 
PYTHIA, UrQMD and HSD events as discussed above. 
In doing so, the acceptance matrices mentioned above are 
used, i.e. the comparison is performed inside the HADES 
acceptance. Furthermore, momenta of leptons are smeared 
in the simulation in order to take into account our finite 
detector resolution. The smearing functions are obtained 
by propagating simulated e + and e~ tracks through the 
detector setup using the Geant package [4*§Il5"0"] , hence tak- 
ing into account the interaction of leptons with the detec- 
tor material as well. This is in particular visible in the sim- 
ulated uj peak shape, where the tail towards low masses 
is due to the energy loss of electrons via electromagnetic 
radiation (bremsstrahlung) . On the other hand, the ion- 
ization (collisional) energy loss of electrons shifts the pole 
position of the reconstructed omega peak by 1 % down- 
wards. 



3.2 Invariant mass distribution 



The comparison of the experimentally measured invariant 
mass distribution of e+e~ pairs to the PYTHIA+PLUTO 



results is presented in Fig. |2(a) The simulated cocktail 
(the sum of the different cocktail contributions is plotted 
as a green curve) reproduces results of GiBUU [25] and de- 
scribes the data reasonably well except for the mass range 
around 0.55 GeV/c 2 where the yield is underestimated. 
The latter deviation is not too surprising, as the virtual- 
ly of the photon 7* — > e + e~ reaches quite high values, 
and therefore it is not guaranteed that the 4(1232) form 
factors fixed at the photon point are still valid. Further- 
more, contributions of higher 4 and N* resonances are 
not includ ed wh ich might also lead to some deficit. 

In Fig. 2(b) we compare our results to the UrQMD [2"2"| 



predictions. In this case, the simulated cocktail clearly un- 
derestimates the contribution from rj and uj mesons, while 
the contribution from the 4(1232) Dalitz and p decays are 
too strong. The latter might be due to the large R— >Np — > 
Ne + e _ couplings as mentioned already above. The small 
rj contribution can be related to the fact that in UrQMD 77 
meson production is mediated by the N*(1535) excitation 
only, which is known to be important for the exclusive 
channels. For the inclusive production, however, channels 




M ee [GeV/c 2 ] 



(a) 




0.2 0.4 0.6 0.8 1 1.2 



M ee [GeV/c 2 ] 



(b) 




M ee [GeV/c 2 ] 



(c) 



Fig. 2. (color online) The HADES data for the p(3.5 GeV) + 
p — > e + e~ X reaction, compared to a simulated cocktail from 
the a) PYTHIA, b) UrQMD and c) HSD event generators. 



6 



HADES collaboration: Inclusive dielectron spectra in p+p collisions at 3.5 GeV 



with one and two additional pions are probably important 
as well. 



Fig. 2(c) illustrates the comparison of the same experi- 
mental data to the HSD results. In this case, the simulated 
cocktail has too strong contributions for all components, 



except the 7r°. The comparison of Figs. 2(a) 2(b) and 2(c) 
points to a lack of understanding of the relative strengths 
of the Z\(1232) and low-mass p contributions. For instance, 
as discussed in 0], implementation of the Iachello model 
|15j enhances dielectron yield in the high mass region. This 
has also been demonstrated in [5T|l25], We expect hence 
that our data can help to clarify the issue of form factors 
in the 4(1232) Dalitz decay and the role of other baryon 
resonances, as mentioned in the introduction. 



3.3 Transverse momentum distributions 

The e + e~ pair transverse momentum p t distributions for 
different invariant mass bins inside the acceptance of HADES 
are presented in Figs. 3(a)[ |3(b)| |3(c) and compared with 
the respective results from PYTHlA+PLUTO, UrQMD 




and HSD calculations. In the mass range M e + e - [GeV/c 2 ] 

< 0.15, dominated by the contribution from dielectrons 
stemming from the ir° Dalitz decay, the experimental p t 
distributions are in reasonable agreement with the simu- 
lated pt distributions from PYTHIA both in shape and ab- 
solute yield. The next mass range of 0.15 < M e + e - [GeV/c 2 ] 

< 0.47 constrains the cross sections of the rj meson and 
the Z\(1232) isobar. In this mass interval, the low-p t part 
of the spectra is populated mainly by the pairs originating 
from ?; Dalitz decays. The high-p t part contains in addition 
a substantial contribution from the 4(1232) Dalitz decay. 
Again, tuned PYTHIA reproduces the experimental data, 
while UrQMD is low by a factor of 2 at low pt and a fac- 
tor of 5 too high at large p t - The HSD results in this mass 
bin overestimate the experimental results by large factors 
at high p t . In the mass interval of 0.47 < M e + e - [GeV/c 2 ] 

< 0.7 the experimental data cannot be satisfactorily de- 
scribed by any of these models. The p t distribution in the 
bin M e + e - [GeV/c 2 ] > 0.7, dominated by direct decays of 
the vector mesons, is again overestimated by the UrQMD 
and HSD event generators. 



3.4 Rapidity distributions 

In contrast to the pt distributions, the rapidity distribu- 
tions do not exhibit such large differences between the 
different models (see Fig. |4| . This may originate from the 
fact that the impact of the large-p t region to these spectra 
is reduced due to the p t integration in each rapidity bin. 



PJGeV/c] 



PJQeV/c] 



0.47 < M^GeV/c 2 ] < 0.7 












y : 




P,[GeV/c] 



PJGeV/c] 



(a) 




PJGeV/c] 



PJGeV/c] 



4 Extraction of the n°, 
sections 



rj, A, p, and tu cross 



(c) 



From the spectrum presented in Fig. 2(a) the inclusive 
cross section of the neutral pion production inside the 



Fig. 3. (color online) Comparison of experimental pt distribu- 
tions to the a) PYTHIA, b) UrQMD and c) HSD events for 
different e + e~ invariant mass ranges as indicated. 



HADES collaboration: Inclusive dielectron spectra in p+p collisions at 3.5 GeV 



7 





M e . e [GeV/c s ]<0.15 


t data 




71 -» e'ey T 






A — } eVN - 






w -> jt°e*e" "I 






P ^ e + e- E 























M E . E [GeV/c 2 ] > 0.7 




0.5 1 1.5 





(b) 





l^. B [GeV/c ! ]<0.15 


t data 




7t° -» e*e"T 1 






A -> eVN " 






» -* n°e + e" 1 






p -i eV 












sum ~= 











(c) 

Fig. 4. (color online) Comparison of experimental rapidity 
distributions to the a) PYTHIA, b) UrQMD and c) HSD events 
for different e + e~ invariant mass ranges as indicated. 



HADES acceptance was estimated from its Dalitz decay 
channel, ir° — > e + e~7, by integrating the yield in the 
mass interval between and 0.15 GeV/c 2 , and by cor- 
recting for the branching ratio (1.198±0.032)% [52]. This 
result is further extrapolated to full phase space using sim- 
ulated events. To investigate the model dependence, the 
extrapolation procedure is repeated for the PYTHIA and 
UrQMD event generators. The obtained inclusive neutral 
pion production cross section in full phase space amounts 
to 16 ± 2.6 mb and 18 ± 2.7 mb, depending on the model 
used for the extrapolation. The errors for this cross sec- 
tion consist of a statistical part (negligible in the pion 
region) and a systematic part stemming mainly from the 
normalization to the elastic p+p cross section and the ef- 
ficiency correction procedure. The cross section for the rj 
mesons were obtained in a similar way by integrating di- 
rectly its co ntribution to the simulated cocktail presented 
in Fig. 2(a) correcting for its Dalitz decay branching ratio 
(7±0.7) x 10~ 3 [S5] and extrapolating to full phase space 
using again UrQMD and PYTHIA+PLUTO simulations. 
Equivalently, one could directly integrate the experimen- 
tal data after subtracting the contributions from other 
sources in the corresponding mass range. In any of these 
cases there is a contribution from the Dalitz decays of 
the Z\(1232) isobar which is also important in the mass 
range relevant for the r\ meson Dalitz decays Q This is 
best seen in the transverse momentum distributions pre- 
Indeed 



scntcd in section 3.3 



3(a) 



as it is seen from Fig 
the pt distribution corresponding to the invariant mass 
range of 0.15 - 0.47 GeV/c 2 essentially fixes the relative 
cont ribu tion of the ij meson and Z\(1232) baryon (see sec- 
tion 3.3 for details). The contribution from the to Dalitz 
decay in this mass range is fixed by its known branching 
ratio and cross section defined by its direct decay channel. 
After getting a satisfactory descriptio n of t he data with 
the PYTHIA event generator (see Fig. 2(a) ), the produc- 
tion cross sections for the vector mesons were obtained 
from their multiplicities in full phase space generated by 
PYTHIA. The acceptance of the spectrometer in the vec- 
tor meson mass range is high in all dimensions of the phase 
space. Therefore, to a large extent, the cross sections for 
vector mesons obtained from PYTHIA are model indepen- 
dent. Our final numbers for the cross sections are listed in 
Table 1. The uncertainties for the cross sections are calcu- 
lated based on experimental errors (statistical + system- 
atic). Fig. [5] shows a compilation of measured production 
cross sections of vector mesons (w, p) in p+p reactions at 
different energies. Similar plots for different pion species 
and the r\ meson are shown in Fig. [6] where we use the av- 
erage values of the cross sections obtained with PYTHIA 
and UrQMD extrapolations listed in Table l(for details 
see the text above). 



Note that both, A + and Z\° states, are taken into account 
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Fig. 5. (color online) Cross sections for the vector mesons u (left panel) and p (right panel), in p+p collisions as a function of 
Mandelstam variable s. Open circles represent inclusive production cross sections, while the full circles correspond to exclusive 
productions [53] ■ The dashed curve refers to the OBE calculations for the exclusive channels [53] . Cross section values obtained 
in this work are depicted in red. 
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Fig. 6. (color online) Cross sections of pions (left panel) and r\ mesons (right panel) in p+p interactions 53,55,56,57,58,59 . 
The dashed lines in case of pions refer to the parametrizations used in HSD transport model, while for the rj meson the line 
corresponds to its exclusive production through the N*(1535) resonance [60]. The full circles in this case illustrate the measured 
exclusive production cross sections of the rj mesons in p+p collisions. Inclusive production cross section values obtained in this 
work are depicted with red open circles. 





7T° 


7] 




P 


u> 


Oi [mb]: 


18± 2.7 (16 ± 2.6) 


1.14+0.2 (0.93+ 0.14) 


7.5 ± 1.3 


0.233 ± 0.06 


0.273 ± 0.07 



Table 1. Inclusive cross sections obtained for different particles. For ir° and rj mesons the uncertainty due to model dependence 
has been included as explained in the text. The values corresponding to the extrapolation with the UrQMD model are presented 
in brackets. 
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Fig. 7. (color online) Invariant mass distribution of e + e~ pairs 
in the r/ meson mass range. The experimental data is fitted 
with a polynomial distribution (magenta curve) by excluding 
the range of m v ±3<r around r\ meson pole mass. The invariant 
mass distribution of e + e~ pairs from simulated direct dielec- 
tron decays of the r\ meson (blue curve) is presented on top of 
the background (dashed-blue curve). 



5 Direct decays of the 77 meson 



Using the data described above we can derive an upper 
limit for the branching ratio of the direct 77 meson decay 
77 — > e + e~. Fig. [7] shows the invariant mass distribution 
of the e + e~ pairs in the 77 meson mass range. The ex- 
perimental data points are fitted with a polynomial back- 
ground function, excluding some range around the 77 me- 
son mass. There is no visible indication of a peak structure 
from the direct decays of 77 mesons. Indeed, performing 
the Kolmogorov-Smirnov test one gets a 90% probabil- 
ity for the consistency of the background function and the 
data points. However, one can still estimate an upper limit 
for this decay process from the data using the method of 
Feldman and Cousins [51]. In order to test the robust- 
ness of the obtained results, five different ranges of the 
background function have been tested. Moreover, the same 
procedure was repeated for different polynomial functions. 
The output of the Feldman Cousins method is an upper 
limit for the signal counts with a 90% confidence level. 
Using this upper limit for the signal counts and the 77 pro- 
duction cross section, reported in the previous section, an 
upper limit for the branching ratio of (4.9+0.7-1.2) x 10 -6 
is obtained. This value is about 6 times lower than the 
most recent value of 2.7 x 10 -5 from [52]. To demonstrate 
how the hypothetical 77 peak would look like, we show in 
Fig. [7] the 77 shape from a simulation on top of the exper- 
imental data, using the branching ratio of 4.9xl0 -6 and 
the production cross section of 1.14 mb from Table 1. The 
significance of the added hypothetical peak amounts to 
5.5 %. 



6 Summary 

In summary, we reported on a dielectron measurement in 
p+p collisions at 3.5 GeV projectile kinetic energy. For 
the first time the inclusive production cross sections for 
neutral pions, 77, oj and p mesons were determined from 
dielectron experimental data. The experimental distribu- 
tions were compared to the results from the PYTHIA, 
UrQMD and HSD event generators, which use different 
physics assumptions to generate the parent hadrons de- 
caying subsequently into e + e _ at this projectile energy. 
With some minor tunes, PYTHIA+PLUTO results de- 
scribe the experimentally observed distributions in a bet- 
ter way than the resonance production picture used in 
UrQMD. We hope that our data stimulate further work 
pinning-down the issue of formfactors in Z\(1232) Dalitz 
decays and shed light on the role of other baryon reso- 
nances in p+p collisions. Using our data, it was demon- 
strated that the upper bound for the direct 77 — > e + e~ 
decay can be improved by a factor of ~6 compared to the 
value quoted in [62] . 
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